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Abstract 

To investigate the viability of the 4th root trick for the staggered fermion determinant in a simpler 
setting, we consider a two taste (flavor) lattice fermion formulation with no taste mixing but with 
exact taste- nonsinglet chiral symmetries analogous to the taste-nonsinglet U{V)a symmetry of 
staggered fermions. M. Creutz's objections to the rooting trick apply just as much in this setting. 
To counter them we show that the formulation has robust would-be zero-modes in topologically 
nontrivial gauge backgrounds, and that these manifest themselves in a viable way in the rooted 
fermion determinant and also in the disconnected piece of the pseudoscalar meson propagator as 
required to solve the U(l) problem. Also, our rooted theory is heuristically seen to be in the right 
universality class for QCD if the same is true for an unrooted mixed fermion action theory. 
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I. INTRODUCTION 



The use of dynamical staggered fermions in lattice QCD simulations has made it possible 



to Obtain reSU .t s wit, oopreceoeoted M gh preclsi oo H 3 B. However, this approach , 
controversial due to the use of the "4th root trick" : A staggered lattice fermion corresponds 
to 4 continuum fermion flavors (nowadays called tastes to distinguish them from the actual 
quark flavors) , so the fermion determinant for each dynamical quark flavor is represented by 
the 4th root of the corresponding staggered fermion determinant. Since this formulation is 
not manifestly a local lattice field theory there is a danger that it might not be in the right 



universality class for QCD. (In fact it has been argued j^] that this lattice theory is necessarily 
non-local but with locality being restored in the continuum limit [5]. 1 ) Because of the high 
stakes, this has become a prominent, hotly debated issue in the lattice community. E.g., it 
has been the topic of 5 plenary talks at the last 4 annual lattice field theory conferences; 
the corresponding proceedings papers , lloL I 111 . Il2j ] can be consulted for reviews from 
various perspectives. 

While the results to date are in excellent, unprecedented agreement with experiment, a 
major question regarding the 4th root trick for staggered fermions is whether it can work 
in situations where chirality is important. This includes in particular producing the large 
mass of the rj' meson where existence of fermionic zero-modes with definite chirality and 
their connection with topological charge of lattice gauge fields via the Index Theorem plays 
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an essential role [3]. 2 M. Creutz has argued against this in a series of papers [17 



18 



12, 



2 if ] based on the fact that the the taste-nonsinglet U(l) chiral symmetry of staggered 
fermions implies properties of the rooted staggered fermion determinant that do not hold for 
a genuine single-flavor fermion determinant. The subsequent rebuttals of these arguments 



22], |23( rely to a large extent on invoking full taste symmetry restoration on the continuum 
limit. However, Creutz challenges whether this can actually occur in a way that correctly 
reproduces nonperturbative effects connected with chirality. In this situation it is desirable 



1 In the free field case (at least for m ^ 0) the rooted staggered formulation corresponds to a local field 
theory already at non- vanishing lattice spacing Q , as was also confirmed numerically [t| ■ 

2 Efforts to calculate the rf mass to high precision with dynamical staggered fermions are currently underway 
[3j ln the meantime, encouraging evidence that this formulation is able to correctly reproduce 



topological aspects of QCD has been given in Ref. 16| where results for the topological susceptibility were 
presented. 
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to have a simpler setting where the same issues arise and where they can be studied more 
explicitly. We provide and study such a setting in the present paper. 

The paper is organized as follows. In sect. II we contrast a general mixed fermion formula- 
tion with a rooted formulation based on a 2-taste lattice Dirac operator without taste- mixing, 
showing heuristically that if the former is in the right universality class for QCD then so 
is the latter. In sect. Ill we introduce the specific 2-taste lattice Dirac operator with exact 
taste-nonsinglet chiral symmetries on which the rooted formulation studied in this paper is 
based. In sect. IV we study the properties of the single-flavor fermion formulations based 
on the 1-taste lattice Dirac operators making up our 2-taste operator, and use this to derive 
properties of the rooted formulation based on the latter. In sect. V we discuss the pseudo- 
scalar propagator in the rooted formulation, and conclude with a discussion in sect. VI. A 
relation between, and differences between, our 2-taste formulation and the 2-flavor Wilson 
fermion theory with twisted mass is discussed in an appendix. 

II. PRELUDE: MIXED FERMION ACTION VERSUS A ROOTED FORMULA- 
TION 

We begin with some general remarks on generating functionals for lattice fermions (in a 
fixed gauge field background, in Euclidean spacetime). For a single quark flavor described 
by a lattice Dirac operator D\ the generating functional is 



For a "mixed fermion action" where the sea quark is described by D\ and the valence quark 
by another lattice Dirac operator D 2 the generating functional becomes 




det D 1 lr > 



(1) 



Z f ,mixed(v,v) = detD 1 e r > ? 



(2) 



Writing this as 



Zf, m ixed(V,V) = e 



det D 2 e^ lr > 



(3) 



where 



AS = Tr log D 1 - Tr log D 2 



(4) 
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we see that the full lattice QCD theory with mixed fermion action is equivalent to the lattice 
fermion being described solely by D 2 and the lattice gauge field action being shifted by 



S gauge * S gauge ~~\~ AS . (5) 

If the shift (JSJ) does not change the universality class, i.e., leaves the lattice theory in the 
right universality class for QCD, then surely the same is true for the smaller shift 

Sgauge * Sgauge ' (^) 

But this shift is equivalent to leaving S gaU ge unchanged and changing the fermion determinant 
in the mixed fermion generating functional (T3J) by 

detA -> (detDxdetA;) 1 / 2 (7) 

We conclude that if the lattice QCD theory with mixed fermion action is in the right uni- 
versality class for QCD then so is the theory where a dynamical quark is described by the 
2nd taste (flavor) of the 2-taste lattice Dirac operator 

CO 



D 2 

and with rooted fermion determinant (det-D) 1 / 2 . 

Normally there would be no reason to consider such a formulation in practice rather then 
just using Di or D 2 or the mixed fermion formulation. But it is useful to consider this 
formulation for theoretical investigation of lattice QCD with rooted fermion determinants - 
it is simpler than the relevant case of staggered fermions since the taste (flavor) interpretation 
is manifest from the beginning and there is no mixing between the different tastes. In the 
next section we will exhibit a 2-taste lattice Dirac operator of the form (|SJ) with properties 
analogous to the staggered Dirac operator and for which Creutz's objections also apply. 
The preceding considerations have already shown (at least heuristically) that the viability 
of using the rooted determinant of such an operator is assured if the related unrooted mixed 
fermion theory is in the right universality class for QCD. 
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III. A 2- TASTE LATTICE DIRAC OPERATOR WITH EXACT CHIRAL SYM- 
METRIES 



The specific 2-taste lattice Dirac operator we will study is given in the massless case by 
D= ( D n n) ' £± = 7„V M ±*7 5 (4A + m 5 ) (9) 



D. 

where a =lattice spacing, V M is the usual symmetrized lattice covariant derivative and A 
the usual lattice Laplace operator. For reasons discussed below we have included a mass 
parameter m§ in the operator. It should not be confused with the usual mass; we introduce 
the latter in the usual way: the massive 2-taste Dirac operator is D + m. Note that D± , 
and hence D, are anti-hermitian, and that in the free field case (link variables set to unity) 

DlD±] f = [yy+ (a±A + m 5 ) 2 l (10) 



free 



free 



which shows that D± is free of fermion doubling so D describes two lattice fermion tastes 
as claimed. Writing D as 

D = ( 7m ® l)V M + i(7 5 ®a 3 )(a|A + m 5 ) (11) 

we see that it has the taste-nonsinglet chiral symmetries 

{D,Vj} = , ^ = 76(8^,^ = 1,2 (12) 

where <7i, a 2 , cr 3 are the Pauli matrices acting on taste space. On the other hand, D 
breaks the symmetry of the chiral transformation generated by 75 <8> 03 (and also 75 <S> 1 
as it should to produce the axial anomaly). Consequently the pion spectrum will not be 
SU(2)-symmetric, so the fermion theory described by D is not equivalent at non-zero lattice 
spacing to a two-flavor theory with both flavors described by the same single-taste lattice 
Dirac operator. 3 

Thus D shares key properties with the massless staggered lattice Dirac operator: it is anti- 
hermitian and has taste-nonsinglet chiral symmetries which protect against additive mass 
renormalization and are expected to be spontaneously broken just like the taste-nonsinglet 
U(l) chiral symmetry of staggered fermions, while other chiral symmetries are explicitly 



3 This is the same reasoning that was used in Ref. [J| to draw the analogous conclusion for the staggered 
fermion theory. 
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broken at non-vanishing lattice spacing. In fact the expression (fTTj) has clear similarities 



with the massless staggered Dirac operator in the flavor (taste) representation 24. |25|. 

However, there is also a significant difference: Our operator breaks the parity and time 
reversal symmetries, since the "Wilson-like" term in D± gives a pseudo-scalar term in the 
fermion action. Consequently, radiative corrections will generate a pseudoscalar mass term. 
Therefore we have included a bare pseudo-scalar mass term with mass in (jUj); it should 
be tuned to a critical negative value as the continuum limit is approached so as to cancel the 
pseudo-scalar mass term generated by radiative corrections and thereby restore the P and T 
symmetries in the continuum limit. (This is analogous to the tuning of the bare scalar mass 
to a critical negative value to reach the chiral limit with usual Wilson fermions.) Through 
this we also expect the chiral transformation generated by 75 £g> cr 3 to become a symmetry 
of the 2-taste theory in the continuum limit. Usually we will suppress the ms-dependence 
of D in the notation, although sometimes we will indicate it explicitly as D(m§). 

The 2-flavor theory described b y D (m?,) + m can be obtained from the 2-flavor Wilson 
fermion theory with twisted mass |26| by a flavored chiral rotation of the fields. We show 
this in Appendix A. However, as shown there, the symmetries of the theories have different 
interpretations: the chiral symmetries ffT2l correspond to vector symmetries in the Wilson 
case. 

A major advantage that our setting has over the staggered one for investigating the 
viability of rooting is that there is a single-flavor fermion theory that our rooted theory can 
be explicitly compared with, namely the the theory described by D + + m (or D_ + m). 
Comparison of the rooted theory based onD + m with the single-flavor theory described 
by D + + m will be our main focus in this paper. Through this we will be able to counter 
Creutz's objections to the rooting trick quite explicitly. 

The starting point for much of Creutz's argumentation against the rooted staggered 
fermion determinant is the observation that, as a consequence of its exact U(l) chiral sym- 
metry, the staggered fermion theory with mass m is equivalent to the one with mass term 
changed by 

mipip — > rrnpip cos(20) + imipT 5 ip sin(20) (13) 

for any 6, where is the operator on staggered fermion fields corresponding to 75 in the 
naive lattice fermion theory from which the staggered theory originates. See Eq.(4) of 
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Ref. [19(1 . In particular, m, —m and ±iT^m are all physically equivalent. Therefore, the 
rooted staggered fermion determinant is invariant under these changes in the mass term, 
unlike the determinant of a genuine one-flavor lattice Dirac operator. Exactly the same is 
true in the present 2-taste theory when T 5 in (fT51) is replaced by our r\ or T 2 in light of (TT21) . 
Therefore, all Creutz's objections against rooting based on (TT5|) apply just as much in our 
case. In the following sections we derive explicit relations between the rooted determinant 
formulation and single-taste theories in the present case which show that, despite Creutz's 
concerns, the rooted formulation does appear to be viable, or at least a good approximation 
when the quark masses are not too small. 

IV. PROPERTIES OF THE SINGLE-TASTE THEORY AND IMPLICATIONS 
FOR THE ROOTED FORMULATION 

In this section we derive general properties of the single-taste theories described by D±+m 
and use them to obtain significant indications of the viability of the rooted determinant 
det(.D + m) 1//2 . Throughout the following we assume that the lattice is finite; then the 
vectorspace of lattice spinor fields is finite-dimensional and the fermion determinants are all 
finite. 

A. Fermion determinants and would-be zero-modes 

D± in (jHD satisfies 

£±75 = -75^t (14) 
implying an equivalence between the eigenvalue equations for D + and D_ : 

D+i/jx = <£> £C(75^a) = -iXMx) (15) 

Using this, we find that the rooted determinant of D + m is given in terms of the eigenvalues 
of D+ by 

det{D + m) 1/2 = l[V\ 2 + m 2 (16) 

A 

Comparing this with the determinant of the single-taste Dirac operator D + 

det(D + + m) =Y[(i\ + m) (17) 

A 
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we see that 

det(£> + m) 1/2 = | det(£>+ +m)| (18) 

This shows that using the rooted determinant is the same as removing the complex phase 
of the single-taste determinant det(-D + + m). We will show further below that the single- 
taste determinant is indeed complex-valued, calculate its complex phase when m is in a 
"chiral region" , and discuss how the complex phase can be removed to arrive at the rooted 
determinant via (TTBl . The would-be zero-modes of the 1-taste and 2-taste theories play an 
important role in this, and we begin by considering them in the following. 

An exact zero-mode would give rise to a factor \m\ in (Tl8l) rather than the factor m which 
would appear in a genuine one-flavor fermion determinant. This difference is expected (see, 



e.g., 
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11]) and is inconsequential as long as considerations are restricted to positive m. 
In practice though we do not expect exact zero-modes for this operator; the most one can 
hope for is approximate, would-be zero-modes that become exact in the continuum limit. 
To produce expected non-perturbative effects it is crucial that there are robust would-be 
zero-modes in topologically non-trivial gauge field backgrounds in accordance with the Index 
Theorem. These are indeed present in this case, as we will now show. 
It is useful to introduce the new gamma-matrices 

% = -«7s7/* (19) 

These form another representation of the Dirac algebra: {7^,7^} = 25^, continue to be 
hermitian (7T = 7^) and have the same chirality matrix as before: 75 = 75. Then D + in (J9j) 
can be expressed as 

D+ = i75(7/,V M + §A + m 5 ) 

= if 5 (D w + m 5 ) (20) 

where is the massless Wilson-Dirac operator constructed with the new gamma-matrices. 
Introducing the hermitian operator 

H(m)=^(D W -m) (21) 

we have 

D+ = iH(-m 5 ) (22) 
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so the solutions to the eigenvalue equation 



H(m)ip\(m) = X(m)tl'\(m) (23) 

give back the eigenvalues and eigenvectors of D + in (IT51) as a special case: A = A(— m 5 ) , 
ipx = ipx{-m 5 ). 

From (I2T1) we see that 

A(m ) = D^^a(^o) = moipxirrio) (24) 

i.e., vanishing of an eigenvalue A(m) at mo corresponds to a real eigenvalue mo of the Wilson- 
Dirac operator with eigenvector ^(mo). It is well-known that the would-be zero- modes of 
the Wilson-Dirac operator are precisely the low-lying real (necessarily positive) eigenvalue 
modes. As (|24|) shows, these correspond to crossings of the origin close to zero (i.e., at 
some small positive value mo) by eigenvalues A(m) of H(m). Furthermore, the low-lying 
real modes of have approximate ± chirality under 75 , and from (12T1) it is clear that the 
sign of the chirality is minus the sign of the slope of A(m) where it crosses the origin at mp . 



This is all well-known and was discussed a long time ago by Itoh, Iwasaki and Yoshie 27]. 
It allows a robust, integer- valued index to be defined for the Wilson-Dirac operator in terms 
of the spectral flow of H(m) in the small m region: it is the difference between the number 
of negative and positive slope eigenvalue crossings. In fact this coincides with the index of 



] . It has been studied numerically in 



27 



291 ]. and analytically 



the overlap Dirac operator 
in 30|, |31[ where it was shown to coincide with the topological charge of the (smooth) lattice 
gauge field in the continuum limit in accordance with the Index Theorem. 

By (I22p - (T2~4")) a real mode of D-^ with eigenvalue m is an exact zero- mode of D + (m^) 
when we set = —mo- For an ensemble of lattice gauge fields generated at sufficiently 
small bare coupling (or with a sufficiently improved lattice gauge action) the low-lying real 
eigenvalues of cluster around a critical (positive) value m c - see, e.g., 27]. We henceforth 



4 The robustness of low-lying real eigenvalue modes and index of the Wilson-Dirac operator is ensured by 
the property that in sufficiently smooth backgrounds the eigenvalues cannot vary arbitrarily under defor- 
mations of the background but are constrained to be close to zero. An upp er bound can be analytically 



derived when the plaquette variables satisfy a bound 1 1 1 — — 1 1 1 < e [3J, [33j . More generally, a bound 
constraining the real eigenvalues to lie in neighborhoods of 0, 2/a, 4/a, 6/a, 8/a can be derived in this 
case [34 1 . 
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tune m 5 so that 

m 5 = —m c . (25) 

Consequently the eigenvalues and eigenvectors of D + are i\{m c ) and i[)\(m c ), respectively. 
Then the would-be chiral zero-modes of are in one-to-one correspondence with would-be 
chiral zero-modes of D + . This is seen as follows. A low-lying real mode of D-^ with eigenvalue 
mo is a zero-mode ip\{ m o) for if (mo) with approximately definite chirality. Since A(mo) = 
and mo is very close to m c it follows that A(m c ) is very small, and the corresponding 
eigenvector ip(m c ) is very close to ip(m ) and therefore has the same approximate chirality. 
Recall that the sign of the chirality is the opposite of the sign of the slope of A(m) at 
m = mo- Therefore, ± chirality corresponds to the sign of A(m c ) being ± if m c < mo and 
=p if m c > m . 

Thus, with m.5 tuned as dictated by (1231) . the low- lying modes of D + are generically would- 
be chiral zero-modes; they are robust since they are tied to the would-be chiral zero-modes 
(i.e., the low-lying real modes) of the Wilson-Dirac operator D^. By (1T5|) the same is true 
for D_. Note that a would-be chiral zero- mode for the Wilson-Dirac operator corresponds to 
would-be chiral zero-modes for each of D + and _D_ with the same chirality; consequently it 
corresponds to two would-be chiral zero-modes for D with the same chirality. Thus we have 
established that D + , D- , and D all have robust would-be chiral zero-modes in sufficiently 
smooth gauge backgrounds, and that the index defined from these equals the topological 
charge Q of the gauge background (or 2Q in the case of the 2-taste operator D) in accordance 
with the Index Theorem, since this holds for the Wilson-Dirac operator. 

The tuning of dictated by ( J25l) is also the appropriate one for restoring in the contin- 
uum limit the P and T symmetries and the chiral symmetry of the 2-taste theory generated 
by 75 <8> CT3- Restoring these symmetries means tuning so that the effective pseudo-scalar 
mass term vanishes. The usual signal for the vanishing of an effective mass term (scalar or 
pseudo-scalar) is divergence of the propagator, and from the discussion leading to (1231) it is 
clear that this tuning produces the desired divergence in the present case. (It is analogous 
to the tuning of the bare mass to a critical negative value to to reach the chiral limit with 
Wilson fermions.) 

Returning now to the fermion determinants, we first note that in approaching the chiral 
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limit the scalar mass m should be tuned such that 



Wow\ « \ m \ « \X non low\ (26) 

where Xi ow refers to the eigenvalues of the would-be zero-modes and X non i w refers to all the 
other eigenvalues. The reason is to achieve appropriate near-chiral limit mass dependence 
in the fermion determinants: from ffl6|) - ffT7|) we see that (126]) is necessary and sufficient to 
get 



det(£>+ + m) « JJ m \{ i\ n oniow (27) 
detOD + m) 1 / 2 « [] H II l^noniowl (28) 

To compare the rooted determinant with the single-taste determinant det(D + + m) we need 
to determine the complex phase of the latter. We now calculate it for m in the chiral region 
fT26|) . starting from 

det(D + +17l) TT 771 TT ■ ^nonlow 



TT '"' TT ^ "nonlow ('29") 

det(D+ + m)| a l m l a |Aw ou ,| 



^low ^-nonlow 



where negligible terms ~ and ~ t— — — have been dropped. To evaluate this we will use 

^nonloiu 

II <=IIH) (30) 

which follows from n a i = 1, a consequence of the fact that the dimension of the vectorspace 
of lattice spinor fields is a multiple of 4. Now recall that the eigenvalues X n0 niow are the 
values at m = m c of X non i ow {m) . Generically these do not cross zero in the small m region; 
in particular they do not cross zero as m varies from to m c . Therefore, 

T~J X non i ow (^Tn c ^j j-j X non i ow (0^j j-j ^/oui(O) (31) 

\ , -^nonZotu (^c) \ , \X n onlow (0) \ , ^iouj(O) 

where the last equality follows from JIa p^jj} = 1, a consequence of the known fact that 
Tr |^|°| = 0; see, e.g., [34]. On the other hand, the eigenvalues Xi ow (m) do cross zero in 
the small positive m region. If Xi ow (0) > then Xi ow {m) has negative crossing slope, corre- 
sponding to a positive chirality would-be zero-mode by our previous discussion. Similarly, 
Xiow(0) < implies a would-be zero- mode with negative chirality. It follows that 

n = ( - ir " (32) 



A, 



OIL' 
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where n± denotes the number of ± chirality would-be zero-modes. This together with (J31 
and (130]) leads to 



[] i " nmiow =r Q (33) 

\ , ^nonlow 

'^nonlow 



where 

Q = n+-n_ (34) 

is the index of the would-be chiral zero-modes and coincides with the topological charge 
in sufficiently smooth gauge backgrounds as discussed earlier. Using this in (I29p we finally 
obtain 



det(L>+ +m) 



The equality becomes exact in the limit — > , t— — > 0, which should be regarded as 

the chiral limit in this setting. The prospects for the possibility of being able to take this 

limit (in principle) are discussed in the concluding section. 

Thus for m in the chiral region (|2"6T) the effect of the complex phase of det(D + + m) is to 

shift the QCD theta- vacuum angle by 9 — > 9 — S | . Since the physical theta- vacuum angle 

I I 

must be zero (or extremely close to zero) [35(, the bare theta- vacuum angle in the lattice 
QCD theory must be chosen such that the shifted one vanishes. This is equivalent to having 
trivial theta vacuum and replacing the fermion determinant (with m > 0) by 

det(£>+ + m) -»■ e^ Q det(L>+ + m) (36) 

which is essentially the same as 

det(£>+ + m) -> | det(D + + m)\ = det(D + m) 1/2 (37) 

in the chiral region. This strongly indicates the viability of using the rooted determinant to 
represent the determinant for a single quark flavor in the present case. 

As a further indication of the viability of the rooted determinant we see from ( 1201) with 
(T55D that 

det D+ = det(ZV - m e ) . (38) 
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Since the Wilson fermion determinant is real and positive, this shows that det(D + m) 1 / 2 
coincides at m = with the Wilson fermion determinant with bare mass tuned to precisely 
the negative critical value which it should have in the chiral limit. Therefore, for very small 
m in the chiral region (126]) . the rooted determinant is very close to the chiral limit of the 
Wilson fermion determinant. 

Flipping the sign of m has the effect of complex conjugation on the phase factor in ( l35l ). 
We note in passing that this is a general property the single-taste fermion determinant: 
From (fITI) and ( fT5l) we easily find 

det(£>+ - m) = det(D + + m)* = det(£>_ + m) (39) 



B. Origin of the complex phase 

For the rooted formulation to be viable, the low energy physics it describes should be 
the same as when the fermion is described by D + + m with a bare theta term included in 
the lattice QCD action to cancel the one produced by det(D + + m). For this to hold, the 
complex phase should originate from the ultra-violet part of the spectrum of D + + m, so 
that it is not a manifestation of low energy aspects of the fermion formulation described by 
D + + m. We show this to be the case in the following. 

In sufficiently smooth gauge backgrounds where the Index Theorem relation between 
chirality of would-be zero-modes and topological charge holds, it is known that for each 



would-be chiral zero-mode of the Wilson-Dirac operator there are 15 "doubler" modes 



341 ] . These are eigenvectors of Dy/ with approximately definite chirality and with large 
(positive) real eigenvalues clustered around specific values: If the approximate chirality of 
the zero-mode is ± then the associated real eigenmodes consist of four eigenvectors with 
eigenvalues « 2/a and chirality =F; six eigenvectors with eigenvalues ^ 4/a and chirality ±; 
four eigenvectors with eigenvalues ~ 6/a and chirality =p, and one eigenvector with eigenvalue 
k 8/a and chirality ± [34]. By (I24|) this implies a corresponding family {'?/'j(m)}j = o i i.... i i5 
of eigenvectors, H(m)ipj(m) = \j(m)ipj(m), with each Xj(m) vanishing at a value rrij close 
to 2p/a for some p G {0, 1,2,3,4}. There are p!( - 4 4 ^, j's for each p, and ipj(mj) has the 
approximate chirality ±(— l) p under 75. 

Noting that H{m c ) can be written for each j as H(m c ) = H{rrij) + (m, — m c )7s we see that 
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each if)j(rrij) is an approximate eigenvector for H(m c ), and therefor also for D + = iH(m c ): 

Him^jirrij) « ±(-l) p (m j - m e )if>j(mj) . (40) 

Since m c « it follows that, generically, the spectrum {iX} of D + contains 15 "doubler" 
eigenvalues associated with each would-be zero-eigenvalue; four of them are ~ ^i2/a; six of 
them are ~ ±z4/a; four of them are ~ =p6/a; and the final one is ~ ±i8/a where ± is the 
chirality of the would-be zero-mode. The contribution of these to the phase factor in (1291) is 

II i^r = (TO W(T<) 4 (±0 = T< (41) 

j=l 

It follows that the total contribution to the phase factor from the doubler modes of all the 
would-be zero-modes is (— i) n+ i n ~ = This reproduces precisely the phase factor in fl33l) . 
which for m > gives the complex phase of the determinant in (I35I) . Thus we have found 
that, at least when m is in the chiral region, the complex phase of det(_D + + m) originates 
entirely from the would-be doubler modes associated with the would-be zero- modes of D + . 
An analogous result holds for det(-D_ + m). The doubler mode eigenvalues are ~ 1/a, so 
we conclude that the complex phases of the determinants are not connected with the low 
energy physics of the fermion tastes described by D + + m and D_ + m. 



C. Determinant phase factor and axial anomaly in the classical continuum limit 

A classical continuum limit version of our determinant phase factor result (|35l) arises as 
a special case of a previous result of Seiler and Stamatescu (SS) |36|. They considered the 
m 5 = case of the lattice Dirac operator 

D e = ^V, + e ie ^(alA + m 5 ). (42) 

which coincides with our D± for 9 = ±7r/2. SS showed that the fermion determinant 
det(De + m) produces a theta- vacuum term e~ l9 ® in the classical continuum limit (with 
m > 0). A simple consequence of their specific result, Eq.(19) of |36(, is 5 

detW + m) eQ 
i-o \det(D£ + m)\ [ 6) 



5 The minus sign in the exponent of e %e ® is erroneously absent in Eq.(19) of 36(; it should be present due 
to minus sign in their Eq.(22). 
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where the gauge background is the lattice transcript of a smooth continuum gauge field 
A (satisfying certain technical conditions) with topological charge Q. For 9 = it/ 2 this is 
obviously a classical continuum limit version of our result (1351) with m > 0. 

The result (I43p is obtained as a straightforward consequence of another result of SS, 
namely that Dg + m reproduces the correct axial anomaly in the classical continuum limit 
for all values of 9. This implies in particular that fermions described by D + + m and 
D- + m both reproduce the correct axial anomaly, so the same is true for the 2-flavor theory 
described by our D + m. We emphasize that both tastes reproduce the correct anomaly with 
the right sign; they do not have opposite signs and so Creutz's concern about cancellation 
of anomalies [19| is not realized here. Although the considerations of SS were without m 5 , 
their results extend almost immediately to m^ ^ 0. This is because m$, just like m, appear 
in the final axial anomaly expression through the dimensionless quantities am^ and am and 
hence drop out in the a — > limit. (However, m plays the role of infrared regulator in 
intermediate stages of the evaluation and must therefore be non- vanishing and positive.) 
Here m*, and m may either be constant or tuned as a function of the lattice spacing as long 
as am 5 — > , am — > for a — > 0. 



V. PSEUDOSCALAR MESON PROPAGATOR IN THE ROOTED FORMULA- 
TION 



We now consider the pseudoscalar meson propagator, more specifically its disconnected 
piece G (x,y), which is supposed to solve the U(l) problem by being non-vanishing in 
the chiral limit in topologically nontrivial gauge field backgrounds and thereby partially 
cancelling the connected piece, resulting in quicker decay, and hence a large mass. This 
cancellation, which was already verified a long time ago in the chiral limit with Wilson 
fermions 27] . requires that G DC '(x,y) in a fixed topologically nontrivial gauge background 
develops a singularity ~ 1/m 2 in the chiral limit, produced by the (would-be) zero-modes 
of the (lattice) Dirac operator. 

For simplicity we restrict to the one flavor case; then, with our 2-taste D we have 

G DC (x,y) = itr[(£> + m)-\x,x)( l5 <g> 1)] ±tr[(D + m)-\y, y){ lh ® 1)] (44) 
We have replaced tr — > 5 tr compared to the usual expression to take account of the two 
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tastes of D. It suffices to consider just one of the factors ~tr[- • •]. A simple calculation using 
(HI- (US]) gives 



|tr[(D + m)- 1 (x,x)( 75 ^l)] =E x2 m 2 ^10*075^0*0 (45) 



Exactly the same expression can be (formally) derived in the continuum irom 0ip\ = ityx 
using the fact that ^(75^) = —iX^ipx). However, in the present lattice setting we do not 
have exact zero-modes in general so (f4"5l) and hence G DC (x, y) vanish at m — 0. (This can 
also be seen directly from the chiral symmetries in (IT2l) since Tj (j = 1,2) commutes with 



7 5 (g) 1. The situation is the same for staggered fermions - see Sect. VIII. F of Ref. 271] .) 
Clearly the m — > limit should not be taken before the continuum limit here. 6 The situation 
is different from Wilson fermions where the chiral limit can be reached by tuning the mass 
to a critical negative value [27]. In the present case, reaching the chiral limit requires being 
able to chose m in the same way as in our discussion of the fermion determinants in the 
previous section, namely, it should be in the "chiral region" (T26l) . Then (T4"5]) becomes 

itr[( J D + m)- 1 (x,x)( 7 5® 1)] « £ - ^ o Jz)7sKn>) ( 46 ) 

which gives the correct chiral limit behavior of ( 1451) and hence also G DC '(x, y). The fact that 
only the would-be zero-modes contribute in (146]) fits well with the observation from previous 
numerical studies that G DC is essentially given by the contribution from low-lying modes 



Wilson-Dirac operator in 



4l| 



this was seen for staggered fermions in 40] and for Wilson fermions using the Hermitian 



From (T46I) and (1281) we see that in the chiral region (T26|) with positive m the weighted 
propagator in the rooted theory, det(-D + ■m) x l 2 G DC (x,y), has the same form and mass 
dependence as obtained from the 't Hooft vertex in the continuum setting. (This is clear, 



e.g., from the description of the latter given in 



21|.) This is clearly not the case for values 



of m which are smaller than specified in (T26l) though. 



VI. CONCLUDING DISCUSSION 



In the rooted fermion formulation based on D + m one would expect that any problem 
connected with chirality would show up most clearly in the "chiral limit" of small bare mass 



The necessity of taking continuum limit before chiral limit is well- known for staggered fermions 37l.l38l.l39| 
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m. We have found no sign of this, having derived quite explicit indications of the viability 
of the rooted formulation when the bare mass is positive and in the "chiral region" (T26l) : 



|^Zoui| ITTi | ^nonlow \ 

(47) 

and in particular in the "chiral limit" 

^ - , — ^ > (48) 

The existence of the chiral region and limit requires a gap in the eigenvalue spectrum of 
D between the eigenvalues {i\[ ow } of the would-be zero-modes and the other eigenvalues 
{i^noniow}- It is plausible that such a gap will open up as the continuum limit (bare coupling 
(7 — > 0) is approached: In this limit the fluctuations of the low-lying real eigenvalues of the 
Wilson-Dirac operator around a critical value m c should become smaller and smaller; then 
the same is true for the fluctuations of {\i ow } around zero when m 5 is tuned to —m c (g) 
(cf. §4). Setting 

fi(g) := max{|A iom |} , f 2 (g) := min{| \ non i ow \} (49) 

we expect 

IM _> o for g - (50) 

h{g) 

Then, tuning m as a function of the bare coupling by, e.g., 

miff) = (A(g)f2(g)) 1/2 (51) 

the chiral limit (|48[) is reached as the continuum limit g — > is taken. This implies that 
a chiral region (j47H exists for sufficiently small (7 (and also at larger g for highly improved 
versions of the lattice actions). 

While the requirement m > for the bare mass has been widely recognized (e.g., in the 
reviews jiol, 



HI]) 7 , we have found here that a more stringent condition is required: 8 

m » max{|A/ ou ,|} (52) 



7 There is however a possibility of extending the rooted formulation with positive m to general complex- 



valued m via the introduction of a theta term, as discussed in the staggered fermion case in Ref. 42 1 



This condition is not at all surprising, and in fact the necessity of it could already be inferred in the 



staggered fermion case from the remarks in sect. VIII. F of 27 1 
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In a lattice formulation of QCD with the fermion determinant for each dynamical quark 
represented by a rooted determinant the dependence of each bare mass m q on g is fixed by 
by renormalization conditions; e.g., by requiring that the lattice QCD theory gives specified 
values for a selection of hadronic mass ratios. In connection with this, Creutz has argued 



43( | that the notion of chiral limit for a single light quark (in practice the up quark) is phys- 



ically meaningless when the other quarks remain massive: He argues that non-perturbative 
instanton effects will produce renormalization scheme-dependent additive corrections to the 
light quark mass. If this is the case then m q = is a scheme-dependent statement for the 
bare mass of the light quark. Then there is no physical reason why the bare mass must 
remain positive in a given scheme (i.e., for a given choice of renormalization conditions) as 
the continuum limit is approached, hence the requirement (|52p may be violated, in which 
case the rooted formulation may fail. On the other hand, if the u and d quarks are taken 
to have degenerate bare mass then the pion spectrum is degenerate and the chiral limit is 
physically well-defined as the limit where the pions become massless. In this case we can 
expect to be able to approach this limit from within the chiral region (14"T1) . This applies not 
only for the present formulation (where the product of the degenerate u and d determinants 
are safely represented by the 2-flavor fermion determinant det(-D + m)) but also for the 
staggered formulation where the determinant product is represented by the square root of 
the staggered fermion determinant. 

Expressions analogous to ( fTTl) for the rooted determinant and ( 1451) for G DC (x,y) hold 
for staggered fermions since the eigenvalues of the massless staggered Dirac operator come 
in pairs The present case is more explicit, since the eigenvalues are those of a 

bona fide single-taste lattice Dirac operator D + , whereas no such origin is known for the 
eigenvalues of the staggered Dirac operator. Nevertheless, the chiral limit issues discussed 
here are the same for staggered fermions. So achieving (14TI) - (|4"8"I) in the staggered fermion 
case is also required for taking the chiral limit there. It is encouraging with regard to this that 
numerical studies with improved staggered fermions find a clear gap in the spectrum between 



the low- lying would-be zero-eigenvalues and the remainder of the spectrum [44], |45|, 146 ]. 

Having seen in §4 that the 2-taste lattice Dirac operator D has robust would-be chiral 
zero-modes in topologically non-trivial gauge backgrounds in accordance with the Index 
Theorem, a natural question is whether the same is true in the case of staggered fermions. 
The numerical studies in 4J, |45] strongly indicate that this is the case. In fact, a version 
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of the techniques used in this paper, supplemented with further calculations, enables the 
robust would-be zero-mode result here to also be established for staggered fermions, thereby 



providing a theoretical basis for the numerical results of 44|, |45[ . This will be presented in 
a forthcoming paper. 

The 2-flavor fermion formulation specified by the lattice operator D introduced here is 
mathematically equivalent to twisted mass Wilson fermions, but the interpretation of the 
symmetries is different: Two of the flavored vector symmetries in the Wilson case correspond 
to the chiral symmetries (fl2|) in our case. Other 2-flavor fermion formulations with flavored 



chiral symmetry have recently appeared 
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inspired by graphene structure. Their 



properties were studied in [49j where a general argument was made that 2-flavor ( "minimally 
doubled") fermion formulations with an exact chiral symmetry must necessarily violate 
parity or time reversal symmetry. The formulation based on D in this paper is another 
example of this: it has two exact (flavored) chiral symmetries and violates P and T symmetry 
due to a pseudoscalar term in the action. 
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APPENDIX A: RELATION TO TWISTED MASS WILSON FERMIONS 

The twisted mass Wilson formulation for two lattice fermion flavors has the action 

Stm = xiltNv + o|A + m + i/i7 5 a 3 )x . (Al) 

Noting that 

i^ 5 (T 3 = -^e-^ 75 ^ , a = n/2 (A2) 
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we see that the flavored chiral rotation of the fields 
leads to 

Stm = ^(jiNfi + z7 5 cr 3 (a§A + m) - fi)ip . (A4) 
This coincides with the action for our 2-flavor theory, 

S = $(D(m 5 ) +m)ip (A5) 

with 

m.5 —> m , m — > —fi (A6) 

At m = our formulation has the exact flavored chiral symmetries generated by Tj = 7 5 <7j 
j = 1,2 (recall ffl2l ). By (lA4j) - (lA6j) m = in our theory corresponds to fi = in the twisted 
mass theory. But from (1A1[) we see that this is just the usual 2-flavor Wilson theory with mass 
m. Thus the symmetries, which in our theory are flavored chiral symmetries, correspond 
in the twisted mass setting to non-chiral symmetries of the usual 2-flavor massive Wilson 
theory with vanishing twisted mass. Specifically, these symmetries, which leave x(Dw + , m)x 
invariant, are the vector symmetries 

for j = 1, 2. 
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